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1. Introduction 

Numerous low-dimensional conductors differ by their layered structure and highly anisotropic elec- 
trical conductivity. Their in-plane conductivity is much higher than that in perpendicular direction, i.e. 
along the OZ-axis. Such an anisotropy is typical of cuprates, e.g. YBaCuO in a nonsuperconducting phase, 
transition metal dichalcogenides (NbSe2, TaS2), graphite and its intercalates in particular, as well as a 
broad family of tetrathiafulvalene salts of (BEDT-TTF)2l3 type (JH]. A common feature of many layered 
conductors is their weak electron-energy dependence on the momentum projection p z . This is due to 
the fact that their interlayer distance a is much higher than the in-plane crystal lattice parameter. Then, 
the wave functions of appropriate electrons weakly overlap, and the energy of single-particle charged 
excitations may be represented by rapidly converging series 

£(p)= E £ «(P-L^ cos ^lT £ ' 

n=0 " 

where the functions e n decrease with increasing n and £i(p_J is much less than co(pj_). Here h is the 
Planck constant, p± = [pn] is the momentum projection onto the layer plane. 

In conductors with high charge carriers density allowance for the electron-electron interaction is 
of great importance. The system of quasiparticles carrying a charge should be regarded as the Fermi 
liquid. According to the Landau-Silin theory the interaction between quasiparticles may be taken 
into account in the form of self-consistent field. In this case, the quasiparticle energy depends on the 
distribution function n{p,r, t) for the other quasiparticles. As a result, the energy of an electron acquires 
the correction 8e to the dispersion law Jl.lt 



/2d 3 v' 
^3«p,p>i(p',r,f), (1.2) 



where ri\ is the correction to the equilibrium Fermi distribution function noie). In quasi-two-dimensio- 
nal conductors, the Landau correlation function L(p,p'), as well as the quasiparticle energy in the "gas" 
approximation tl.l) , depends weakly on p z . 
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The specificity of quasi-two-dimensional electron energy spectrum provides the most favorable con- 
ditions for a detailed study of the electron-electron Fermi liquid correlations. If the initial dispersion 
law dl.lt is unknown, investigation of nonequilibrium processes in the electron system in a stationary 
and uniform external field does not provide novel information about electron-electron Fermi liquid cor- 
relations. The solution of the inverse problem of restoring the electron energy spectrum from experi- 
mental data makes it possible to determine the dependence of the electron energy e(p) = £(p) + Se[p) 
renormalized by the Fermi liquid interaction, at the Fermi surface only. 

However, wave processes are very sensitive to the form and magnitude of the Fermi liquid correla- 
tions. The gas approximation does not provide for propagation of the spin waves in nonmagnetic metals, 
predicted by Silin |8], as well as for a series of Fermi liquid effects. Spin waves and transverse electro- 
magnetic waves in the Fermi liquid have been studied in detail not only in metals (see, for example, the 
monograph OIES) but also in low-dimensional conductors fllUl7h . 

Here we consider longitudinal oscillations in the quasi-two-dimensional Fermi liquid. Following Vla- 
sov flill and Landau flill we should solve the set of equations including the Poisson equation 

divE = 47rp', (1.3) 

the continuity equation 

dp' 

— + divj = (1.4) 
dt 

and the constitutive equation relating the electric field density J to the electric field E ~ exp(-i&)f) with 
allowance for the Fermi liquid correlations. 

For the sake of brevity in calculations, we shall confine ourselves to the simplest case where the elec- 
tron dispersion law contains only zero and first harmonics and the expression dl.lt takes the following 
form: 

£(p) = £o(P_l)+ fj_cos0. (1.5) 

Here £o(p_i_) is an arbitrary function, t±_ is the constant value, 8 = apjl~ x . 

We consider the case where the electric field and the wave vector k are parallel to the normal to 
the layers, i.e. to the hard direction for the electric current. We shall derive the dielectric function with 
allowance for the Fermi liquid interaction (FLI), ascertain the dispersion law for longitudinal plasma 
oscillations and find out the distribution of the high-frequency electric field in a half-infinite sample. We 
shall also consider low-frequency collective excitations in a layered conductor which are possible in the 
presence of two bands of the type (1.5) in its electron spectrum. 



2. Dielectric function 



In order to determine the relation between the current density 

) 



2d s p de 

ztt e — n{p,r, t) 

(27ih) 3 dp F 



and the electric field of the wave, it is necessary to solve the kinetic equation 

W coll {n}. 



dn dn dr dn dp 
dt dr dt dp dt 



Naturally, the collision operator Wc n{rc} vanishes working on the equilibrium Fermi function 

-l 



n (E) ■■ 



1 + exp 



T 



(2.1) 



(2.2) 



(2.3) 



which depends on the energy £ = e(p) + 8e{p,r, t) renormalized by the Fermi liquid correlations. In the 
case of a weak perturbation of the system of charge carriers, i.e. when the value of the electric field of 
the wave is small, the collision integral is the linear integral operator working on the function 



n\(p,r, t) = n(p,r, t) - no[e) = -3>(p,r, t) 



dn {£) 
de 



(2.4) 



13704-2 



High-frequency longitudinal oscillations 



We shall confine ourselves to the T-approximation for the collision integral, where W co \\ corresponds 
to the operator of multiplication nonequilibrium part of the distribution function for conduction elec- 
trons by the the frequency of their collisions 1/t: 

W coll {n] = - 0(p,r, t) . (2.5) 

t oe 

In the kinetic equation, it should be taken into account that besides the electric field, conduction 
electrons experience the force of the self-consistent field of interacting quasiparticles 

do d 

— =- = F = eE- — <5£(p,r, t). (2.6) 
df or 

The current density depends only on the nonequilibrium part of the distribution function of charge 
carriers 

Ur,f) = -e f m*de [ dB [ dw ""^ vOfofl.y.r, t) — (2.7) 
J J J oe 2n i n^a 

Here m* is the cyclotron effective mass, q> is the variable of integration over charge carriers states with 
constant values of energy and momentum projection on the normal to the layers, v= de(p)/dp. 
The kinetic equation for the case under consideration takes the form 

kv z ®-ti) y P-hieEv z = i— , (2.8) 

T 

where *P = fl> - 8e, the relaxation time r is supposed to be sufficiently large. 

Within the chosen model for energy spectrum dl.St . electron velocity across the layers v z = v\ sin(0), 
with v\ = -atxhT 1 , is dependent on 8 only and much less than the characteristic Fermi velocity for the 
motion along the layers. Then, in order to find out the dielectric permeability and finally the spectrum of 
longitudinal vibrations, the averaging over cp of the nonequilibrium distribution function of the charge 
carriers along with the Lanadau-Silin correlation function is sufficient. 

In accordance with the symmetry of the problem, the Landau correlation function L(p,p'), which 
connects the effective and the real distribution of charge carriers, may be represented in the form 

L(0,0') = L o + 2LiCOs(0-0'). (2.9) 

In this case 

n 

/dO F 
— w.e'ymo') (2.io) 
2n 



and one can easily solve the equations <2.8MPuV In particular, for = < V F)/<1>, averaging over 6 we 
have 

ieE wW 

\jj f 9 1 1 1 

°~ k a>+ ^[iT-i+w^o + Aw 2 /^ 2 !;^] ' 



where 



i kv z \ , & 
Wik,co) = (- ^ (lr'^l- . (2.12) 

,kv z + 0)/ l^2_ k 2 v 2 



2L\ i 2 [ dno 

, &> = &>+-, <1) = — -r- 

(l + 2Li) t {2nhf 



f dn 3 



With regard to H2.7t and H2.8t . the dielectric function e = 1 + AnQV)likE takes the form 



e{k,a>) = 1 + 



K 2 (bW 



k 2 

(D+W 



- Lq+Aoj 2 



(2.13) 



where 

K 2 ^Ane 2 (\)^-^- 
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is the square of the decrement of the static screening in the gas approximation. As is easily seen from 
1 12.13) and the dispersion equation e[k, w) — 0, the decrement with allowance for FLI equals k ( \/l + Lq) 1 . 
For the activation frequency of the plasma oscillations we have 

„ 2me 2 v 2 4nNe 2 
0)1= r 7 i (l + Li) = n U + Li), (2.14) 

where 77 = {v\l v-p) 2 . Due to the smallness of 77, the plasma frequency in the quasi-two-dimensional con- 
ductor essentially reduces in comparison with the case of a conventional metal. Making use of the for- 
mula <2.14> we obtain the dispersion law for longitudinal plasmons: 



a) 2 (k) = -4- 



X 

(l + Io)* 2 



fc2 - - "-- iLWrr* 



2k2 k+ xk + ^k-- 2] + XK 



(2.15) 



K(k) = 1 + 

K" 

When the Landau correlation function has only a zero harmonic, the dispersion relation becomes notice- 
ably simpler: 

ev\ /2m k 2 + {\ + Ln)k 2 

0){k) = — -\ =. (2.16) 

n \ a K v /jc 2 + (l/2 + L )fc 2 

3. Distribution of the electric field in a specimen 

Having calculated the dielectric function for the infinite specimen, we can solve the boundary prob- 
lem for the half-space z»0. Landau 119] was the first to consider the penetration of the high-frequency 
longitudinal electric field into the Maxwell plasma. Propagation of the longitudinal waves in a degener- 
ated isotropic conductor was studied in [20!l where a detailed analysis of the boundary conditions for an 
electron distribution function was made. In the layered conductor, the angles of incidence of electrons 
on the boundary do not exceed v\ I « 1, which allows us to suppose that charge carriers are reflected 
specularly. The electric field outside the conductor (i.e. between the capacitor plates, one of which is the 
specimen) is supposed to be given by E{z *S 0) = (0,0, £b exp(-iwf)). Then, the field distribution in the 
specimen is described by the following expression 

00 

£ f dfc 

= — — — exp(ifcz). (3.1) 
m J ke{k,w) 

—00 

The k = pole residue gives the field value in the bulk of the conductor E{oo) which is much less than Eq, 

Eg (*>l 

£(oo) = — , e(0,cu) = l ■ . (3.2) 

e[k,a>) a>{a> + 1Z4 /t) 



if a) « Wp : 



However, at iot » 1, the field attenuates in a nonmonotonous manner because the integral <3.1> contains 
an oscillating component 



4EqF 2 f dxx V x 2 - 1 exp [ik\ xz) 

= — — — : — , r-^r (3.3) 

71 J (x 2 -l 



( [x 2 - 1) [2 + (1 + L )F 2 x 2 Y + (2 + L F 2 x 2 Y 



originating from the branching point k\ = S)lv\. At z » V\I\G)\, the ratio of the oscillating component to 
E{oo) takes the form 

Eliz) l-F 2 (Vi\3l2 



— exp(izw/vi). (3.4) 
V S>z I 



Eoo V2n(l + L F 2 /2) 2 [ cbz> 
Here F = ti)/a) p . In the layer of the order of the electron free path l z = v\T, the field oscillates with the 
period 2nv\lw. In the layered conductor, both scale lengths reduce in comparison with the isotropic 
metal, remaining, however, microscopic for non-high frequencies w and t -1 . 
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For the sake of completeness, the resonance case to = co p should be considered. At such high frequen- 
cies, the direct stimulation of the monochromatic field is hard to achieve, but the character of the spatial 
distribution of the resonance harmonics may become apparent in the pulsed mode, as well as in the 
experiments with electron beams. 

Near the plasma frequency co p , the radicals in the expression for the dielectric function may be ex- 
panded into a series about small k 2 , and in the main approximation we have 

(til 

e(0,«) = l , |e(0,w)|«l. (3.5) 

LO{tO + iL\lj) 



4. The spectrum of the short-wave plasmons 

Consider the case of large values of k in collisionless limit. Strictly speaking, this is the case where the 
quasiclassical formula for the conductivity 

o-(k,w) = -ie 2 / — ^— \ (4.1) 
\ kv- a) I 

is inappropriate because it is not invariant with respect to translation k —■ k + 2nla. The translation- 
invariant dielectric function e{k) - l + 4ma{k)/to can be easily constructed if we note that the magnitude 
Tzkv is the result of the expansion in small k of the energy difference in the quantum perturbation theory 
formulas. In the case of the spectrum <1.5> , this difference is 



e L +n _ \ e L n _ \ =hv I U(t|) 



i.e., in (4.1) k should be replaced by Q = {21 a) sm{kal2). As a result, we have 



K 2 



e{k,ct>) = 1 H - 

Q2 



( \ 

(!) 

1 ; 

io 2 -Q 2 v\ 



(4.2) 



Usually, the application of the quasiclassical approximation for the calculation of the plasmon spectrum 
is limited by the condition hk « p^,h/a (see, for example, 1 21]). In the case under consideration, this 
restriction is of no necessity because the dispersion law <1.5> is determined everywhere in the Brillouin 
zone. From <4.2> we obtain: 

[2 b+sin 2 {ka/2) 
to(k) = to p0 \ - — (4.3) 

v b V 2b + sin 2 (ka/ 2) 



where w p o = {evilH) \/2mla, b = {ko) 2 /4 = ane 2 lh 2 i.e., the ratio of the lattice period to the Bohr radius. 
As far as Bohr radius (at m = mo) is of the order of 10~ 8 cm, the parameter b may be very large, especially 
in artificial superlattices with great separation between conducting layers. In this case, the spectrum (473} 
corresponds to a narrow band 



CO 



max 
P 



b+l f 3/4b (b» 1) 

1= , . PTTT ■„ ' (4.4) 



too Vb 2 + b/2 I V2/b {b«l) 

and the static screening decrement k is replaced by 

k = — arsinh ( . (4.5) 

It is easily seen that the results obtained above are valid in an external magnetic field applied across 
the layers, because it has no effect on the electron velocity along the normal to the layers for the energy 
spectrum under discussion. 
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5. Two-band model 



The Fermi surface is supposed to be singly connected. However, the electron structure of "synthetic 
metals" is quite complicated. It is probable that there exist two charge carrier groups (for example, elec- 
trons and holes) of the type dl.St with the different values V zmax : v[ 2) = Bv^ and B 3= 1. In the case of 
two charge carrier groups in an isotropic metal FLI (even in its simplest form L(p,p') = -L a ,/3, a,fi,= 1,2), 
there appears a longitudinal collective mode of the type of zero sound at sufficiently low frequencies, 
t -1 « to « nip. Its phase velocity V = uv^ satisfies the characteristic equation 



D{k) = J^( a W a = LW 1 W 2 = (iW(.Bu)+( 2 W(u)+LW(Bu)W(.u)=0, 



(5.1) 



where 



Ca- 



in,, 



L - ClC2 (^11 + £22 - - L2l) 



m\ + 1112 

and W{u) = W(a)/kvf^) is given by the formula (2.12) . The equation (5.1) has a real root {u > 0) only at 
sufficiently intensive FLI, namely at 



L>L wia iB)=-( 2 /WiB). 

In the quasi-two-dimensional metal, the threshold value L m ; n is much less than that in the isotropic metal 
for the same B, especially near B = 1. In this case L m [ n = 0, and the solution of the equation (5. It is 



u(L) 



l + L 
VI+2L' 



(5.2) 



The wave under consideration corresponds to antiphase partial oscillations of the electron densities 
for both charge carrier groups. The elastic force generated by non-equilibrium charge carrier distribu- 
tion allows one to detect the electron zero sound by means of the concomitant elastic wave. This effect 
was observed in a series of ordinary metals (W, Al, Ga) 12111 . It may be described by jointly solving the 
electron kinetic equation and the elastic theory equation for the half-space with the given oscillation of 
the boundary woexp(-iwf). As a result, for the displacement of ions in the ^-representation, we obtain 



um 

Rik) 



1 + 



2U 
ik 

wA 2 {l) 

(bps 2 



s 2 k 2 {l + R)-co 2 



(5.3) 



1- WiW 2 



l + L 

mk) 



where p is the mass density, s is the velocity of the longitudinal sound. We assume that s « V\. In this 
case, U{k) has two different poles: the pole k s = tol s is connected to the ordinary sound; the other one, 
k\ s , is close to the root of the equation J5.lt and describes an extra elastic wave. It is easy to see that in 
the degenerated case (l/J 1 ' = vf' = v) we have: 



wA 2 (l> 21 
Uos - U — j- (1(2 - — - exp 



copv 2 



l + L 



dzVl + 2L 
Vid + L) 



(5.4) 



This expression is proportional to v l 2 instead of u F 2 in an ordinary metal. 



6. Conclusions 

The analysis given above shows that the electrodynamic characteristics of the layered conductor es- 
sentially differ from those of the isotropic metal with the same charge carrier density. In particular, the 
activation frequency and the velocity of the longitudinal waves, propagating along the normal to the lay- 
ers, decrease essentially. Moreover, the case of quasi-two-dimensional charge carrier spectrum is most 
favorable in observing the electron zero sound and in studying the electron phenomena in various mod- 
ifications of carbon. 



13704-6 



High-frequency longitudinal oscillations 



References 

1. Wosnitza J., Fermi Surfaces of Low-Dimensional Organic Metals and Superconductors. In: Springer Tracts in 
Modern Physics, Springer, 1996. 

2. Singelton J., Rep. Prog. Phys., 2000, 63, 1111; doi 10.1088/0034-4885/63/8/201 

3. Kartsovnik M.V., Chem. Rev., 2004, 104, 5737; doi 10.1021/cr0306891 

4. Kartsovnik M.V., Peschansky V.G., Fiz. Nizk. Temp., 2005, 31, 249 [Low Temp. Phys., 2005, 31, 185-202; 
doi [TuTl063/l . 18844221 . 

5. The Physics of Organic Superconductors and Conductors, ed. by A.G. Lebed, Springer Series in Material Sciences, 
Springer Verlag, Berlin, 2008. 

6. Landau L.D., Zh. Eksp. Teor. Fiz., 1957, 32, 59 [Sov. Phys. JETP, 1957, 5, 40]. 

7. Silin V.P., Zh. Eksp. Teor. Fiz., 1957, 33, 495 [Sov. Phys. JETP, 1957, 6, 387]. 

8. Silin V.P., Zh. Eksp. Teor. Fiz., 1958, 35, 1243 [Sov. Phys. JETP, 1958, 8, 871]. 

9. Pines D., Nozieres Ph., The Theory of Quantum Liquids, Vol. 1, Benjamin, New York, 1966. 

10. Abrikosov A.A., Fundamental Principles of the Theory of Metals, Vol. 1, New York, 1988. 

11. Gokhfeld V.M., Peschansky V.G., Fiz. Nizk. Temp., 1999, 25, 43 [Low Temp. Phys., 1999, 25, 32; 
doi 10.1063/1.593703]. 

12. Kirichenko O.V., Peschansky V.G., Fiz. Nizk. Temp., 2001, 27, 519 [Low Temp. Phys., 2001, 27, 380; 
doi 10.1063/1.13747241. 

13. Peschansky V.G., Stepanenko D.I., Phys. Met. Metall., 2001, 92, 119. 

14. Peschansky V.G., Stepanenko D.I., Pisma v Zh. Eksp. Teor. Fiz., 2003, 78, 772 [JETP Lett., 2003, 78, 322; 
doi [mil34/l. 1 625734]. 

15. Kirichenko O.V., Peschansky V.G., Stepanenko D.I., Zh. Eksp. Teor. Fiz., 126, 1435 (2004). [J. Exp. Theor. Phys., 
2004, 99, 1253; doi |10.1134/lTl 854813 1. 

16. Kirichenko O.V., Peschansky V.G., Stepanenko D.I., Phys. Rev. B, 2005, 71, 045304; 
doi 10.1103/PhysRevB.71.045304 

17. Kirichenko O.V., Peschansky V.G., Stepanenko D.I., Condens. Matter Phys., 2005, 8, 835. 

18. Vlasov A.A., Zh. Eksp. Teor. Fiz., 1938, 8, 201. 

19. Landau L.D., Zh. Eksp. Teor. Fiz., 1946, 16, 574. 

20. Gokhfeld V.M., Kaganov M.I., Lyubarskii G.Ya., Zh. Eksp. Teor. Fiz., 1987, 92, 523 [Sov. Phys. JETP, 1987, 65, 423]. 

21. Bezugly E.V, Burma N.G., Deineka E.Yu., FiT V.D., Sverkhprovodimost', 1991, 4, 661 (in Russian). 

BMCOKOMaCTOTHi n03flOB>KHi OCLJI/l/lflLjiT KBa3i-fl,BOBMMipHOl' 

esiexTpoHHOi pifli/mi/i 

B.M. rox^e/ib/P 3 , O.B. Ki/ipi/meHKcP, BJ. niu^aHCbKnP 

1 /],OHeL4bKHti (|)i3HKO-TeXHiHHwR iHCTHTyT iM. A.A. Ta/lKil-ia, flOHei4bK, YKpaiHa 

2 ct)i3WKO-TexHi4Hnti iHCTHTyT HH3bKnx TeMnepaTyp iM. B.I. BepKiHa HAH YKpai'HH, XapKiB, YKpaiHa 

npoaHa/ii30BaHO cneu,w(f>iKy no3flOB>KHix ko/ickthbhux e/ieKTpoMamiTHHX KO/inBaHb y LuapyBaTOMy npoBiflHW- 
Ky3 KBa3i-flBOBHMipHHM e/ieKTpOHHUM eHepreTUHHUM cneKTpoM. 
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